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1. Features 

We present ELSA, a new modular software package, written in C, to analyze and man- 
age spectroscopic data from emission-line objects. In addition to calculating plasma diagnos- 
tics and abundances from nebular emission lines, the software provides a number of conve- 
nient features including the ability to ingest logs produced by IRAF's splot task, to semi- 
automatically merge spectra in different wavelength ranges, and to automatically generate var- 
ious data tables in machine-readable or I^TJ^X format. ELSA features a highly sophisticated 
interstellar reddening correction scheme that takes into account temperature and density ef- 
fects as well as He II contamination of the hydrogen Balmer lines. Abundance calculations are 
performed using a 5-level atom approximation with recent atomic data, based on R. Henry's 
ABUN program. Downloading and detailed documentation for all aspects of ELSA are available 
at ,http: / / www.williams.edu / Astronomy / research/PN / ^ 

1.1. Input 

ELSA as a tool for astrophysicists has a two-fold purpose: it automates some of the more 
laborious aspects of spectral data analysis for emission-line objects, while also providing highly 
accurate and refined facilities for dereddening corrections and atomic abundance calculations. 
Furthermore, it can generate output to screen, to spreadsheet delimited file formats, or to 
M5h;X source files for each step of the analysis process. 

ELSA takes the logfile output from IRAF's splot task and matches the measured wavelengths 
(adjustable for a user-input radial velocity) to a list of canonical reference wavelengths provided 
by the user. ELSA can then perform a variety of tasks, including interstellar reddening correction, 
generation of a wavelength-flux data table, computation of plasma diagnostics and elemental 
abundances, and flagging unknown or commented wavelengths for later review. User-deflned 
comments in the original logfiles are preserved throughout the analysis process. It is also possible 
to bypass the reddening correction phase in order to run analysis of data that has been previously 
dereddened independent of ELSA. 

1.2. Reddening Correction & Tg-Nf, Loops 

Dereddening is performed through an iterative process that works in several stages: correcting for 
the He^^ Pickering line contamination of the Ha and H/3 lines, accounting for the co-dependence 
of Te and Ne determinations, dynamically calculating the appropriate Ha/H/3 ratio based on 
Te and Ne, and calculating the reddening coefficient c based on the appropriate intrinsic Balmer 
ratio. ELSA makes a first-pass correction assuming a user-defined default Hq/H/3 {e.g., 2.86). 
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This rough value is used to apply a reddening correction from Savage & Mathis (1979) to Ha 
and to He II A4686. The corrected A4686 flux is then used to calculate and then subtract the 
portion of the fluxes at Ha and H/3 that come from He"'"'^ Pickering line contamination. The 
reddening correction is then run inversely to "reredden" Ha and H/3. A new value of c is then 
computed and the process repeats until successive values of c converge. The final c value is used 
to deredden [O III] A4363 and A5007, which are used to derive Te assuming a fixed "seed" Ne . A 
2D polynomial interpolation table is used to determine Ne using the previously calculated Te. 
A new Te is then calculated based on the improved Ne, and so on until both approach stable 
values. 

The final values for Te and Ne arc used to refine the intrinsic Ha/H/3 ratio, using data and 
code based on Storey & Hummer (1995). ELSA then returns to the beginning of the loop and 
malces a second pass at dereddening the raw data with the new Ha/H/3 used to determine c. 
The entire loop continues to iterate in this manner until Ha/H/3 itself converges satisfactorily. 
The final c value is then used to deredden all lines read from the raw data. ELSA then stores 
this as the pool of data to be used for generating tables or computing chemical abundances. 

1.3. Chemical Abundance Calculations 

ELSA uses a standard five- level atom model as described in Kwitter & Henry (2001), with 
code based on R. B. C. Henry's ABUN routine. ELSA can compute abundances for various 
ionization stages of He, N, O, S, Ar and Ne; total abundances are computed using ionization 
correction factors as described in Kwitter & Henry (2001). Atomic data are drawn principally 
from Osterbrock & Ferland (2005) or from Mendoza (1983). Currently, ELSA contains a two- 
zone ionization model, with T[i^ jj] and T[o j/j] regions. 

1.4. Tabular Output Options 

ELSA can produce output from its tasks in variety of formats. The dereddencd line fluxes can 
be output to a tab-delimited file, a comma-separated value (CSV) file readable by spreadsheet 
programs such as Microsoft Excel, or a MJgX source file suitable for insertion into a publication 
manuscript. The content of this table is controlled by a user-defined list of wavelengths, and 
corrected line intensities are given relative to H/3=100. ELSA also has the capability to produce 
a composite table of the same lines across multiple objects in M^jX format, a product that 
would otherwise be quite laborious to produce. Comments read from log files are also preserved 
and output into the tabular results 

1.5. Planned Improvements 

Future plans for ELSA include the incorporation of robust error propagation, the addition of 
UV and IR analyses, expansion to a multiple-zone model (similar to IRAF's nebular package) 
and support to accommodate extragalactic emission-line objects. 
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